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Abstract 

We apply the canonical quantum field theory based on the Liouville-von Neumann equation to 
the nonequilibrium linear sigma model. Particular emphasis is put on the mechanism for domain 
growth of disoriented chiral condensates due to long wavelength modes and its scaling behavior. 
Scattering effects, decoherence and emergence of order parameter are also discussed beyond the 
Hartree approximation. 
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I. INTRODUCTION 



The high density and temperature state of hadronic matter consists of quark-gluon plasma 
and would have occurred in the early universe or may be realizable in heavy ion collision 
experiments. In QCD with two massless quarks, the chiral symmetry SU{2)l x SU{2)ii 
at high temperatures spontaneously breaks down to SU(2)l+r at lower temperatures by 
the quark- ant iquark condensate {q\qRj) = crSj + ivr ■ 7^, an order parameter. The field 
(pa — (c, tt) respects 0(4) rotations and thus belongs to the universality class of a four 
component isotropic Heisenberg antiferromagnet [1]. The effective theory of the QCD phase 
transition is described by the linear sigma model [2] or equally by a mean-field theory of 
Polyakov loops and/or by the glueball fields for the hadronic states of QCD [3]. 

In QCD the SU (2) phase transition would probably be a second order whereas the SU (3) 
phase transition would be a weakly first order. In relativistic heavy ion collisions, hot and 
high dense regions are made, where the chiral symmetry would be restored, and as the regions 
cool, the quark-gluon plasma would undergo the second order phase transition. In second 
order phase transitions, as temperature approaches a critical temperature, the correlation 
length cannot grow indefinitely and must be frozen due to causality [4]. However, the 
cooling process may be fast enough for the kinematic time scale of quench to be smaller 
than thermal relaxation time scale. The rapid expansion of the quark-gluon plasma enforces 
a rapid quench and results in a phase transition far from equilibrium (nonequilibrium) . In 
such a nonequilibrium phase transition domains (regions of misaligned vacuum) develop due 
to the instability of long wavelength modes. 

Another possible candidate for the QCD phase transition is Centauri events in comic 
rays [5]. Anomalously large event-by-event fiuctuations have been observed in the ratio of 
charged to neutral pions and thus require a new theory or interpretation. The disoriented 
chiral condensate (DCC) of classical pion fields was introduced as one of the proposed 
mechanisms for coherent emission of pions from a large domain [6]. If the QCD phase 
transition proceeds in equilibrium, all directions of tt are equally probable and domains of 
size 1/Tc cannot explain the anomalous production of pions of a certain kind. Rajagopal 
and Wilczek advocated the nonequilibrium QCD phase transition through a quench for DCC 
domain growth, where unstable long wavelength modes of pions are exponentially amplified 
[2]. 
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In this talk, we adopt a recently introduced canonical field method to elaborate the 
nonequilibrium QCD phase transition based on the linear sigma model. The quench process 
is imitated by a mass squared that changes signs during a finite quench time [7, 8, 9]. 
We particularly focus on the dynamical process of domain formation from long wavelength 
modes of pions growing exponentially during an instability period and on the scaling behavior 
of domains size. The scaling behavior of domains has been found through simulations in 
condensed matter systems and cosmology. We further discuss the effects of direct scatterings 
among each pion field modes on domains size and on decohering long wavelength modes and 
emergence of an order parameter. 



II. NONEQUILIBRIUM LINEAR SIGMA MODEL 

The effective theory of QCD phase transition with two massless quarks is described by 
the quark-antiquark condensate {qq), which, in turn, defines (pa — (c, 7?), 7? being the pion 
field. The quark mass niq provides a symmetry breaking external field to an, otherwise, 
(9(4) symmetric field theory at high temperatures. To include a quench, we may write the 
linear sigma model [10] in the form 

1. , A. ,„ , 1 



(1) 



The time- dependent quench process has been imitated by the symmetry breaking mass 
squared m^it), which changes signs during the phase transition. In this sense Eq. (1) is 
the nonequilibrium linear sigma model. As H does not affect much the instability of long 
wavelength modes during the phase transition, we assume if = for simplicity reason. 

To illustrate how nonequilibrium phase transitions in general affect formation of domains, 
we consider a simple field model with a quench time scale [7]. The simple model motivated 
by the nonequilibrium linear sigma model is given by the potential 

V{4>a)^'^r4>a+\{r4>a)\ (2) 



with the mass squared 



m^{t) = m\ — rngtanh^^— ^, (mg > m^). (3) 



In the limit of zero quench time (r = 0), we obtain the instantaneous (sudden) quench. 
Now finite temperature field theory [11] cannot be applied to this model when the quench 



time scale r is smaller than relaxation time scale. Further, in the second order phase 
transitions, long wavelength modes grow exponentially during rolling over the barrier. Finite 
temperature field theory does not properly take into account dynamical processes of phase 
transitions. 

One therefore needs some nonequilibrium quantum field theory when the kinematic time 
scale is smaller than thermal relaxation time scale so that finite temperature field theory 
cannot be applied. There are several methods for nonequilibrium quantum fields such as the 
closed-time path integral [12] and the functional Schrodinger-picture [13]. Recently there has 
been developed a canonical method based on the quantum Liouville-von Neumann equation, 
which provides all time-dependent Fock states at the leading order [7, 8]. The new canonical 
method is equivalent, at leading order, to the time-dependent Hartrcc approximation, but 
it can go beyond the Hartree approximation since any perturbation method can be readily 
applied to these Fock states. 

As a nonlinear theory, the linear sigma model has defied yet any nonperturbative solution 
in a closed form. We can, at best, rely on perturbation methods. The Hartree approximation, 
though being a perturbation scheme, includes some part of nonperturbative effects at the 
lowest order [14, 15]. In the Hartree approximation, dividing the a field into a background 
field (f){t) and its quantum fiuctuation x(x, i), and using the Hartree factorization [14], we 
obtain the truncated Hamiltonian for the linear sigma model 

/ 4 ^ ^ + ^ ^ + + !^.-] , (4, 

where the effective couplings are 

h^{t) = (/.(i)[m2(i)+4A02(t)+4A(7f2)], 

mlit) = m2(t)+4A02(t)+4A(7f2), 

ml{t) = m\t) + 12X(f>\t) + 4A(7r2). (5) 

III. DCC DOMAIN GROWTH 

Quantum dynamics of DCC can be further approximated by an exactly solvable model 
now motivated by the linear sigma model (1) or the Hartree approximation (4). The model 
Hamiltonian for the pion field (f)a is [16] 

Ho{t) = J d'x[^n''na + ^{Vr -^(Paf + lm''{t)r<l>a]. (6) 



Note that all pion fields are decoupled from each other since the nonlinear term is neglected 
at this moment. In terms of the Fourier cosine and sine modes 

and with a compact notations a — {(±), k}, the Hamiltonian becomes a sum of decoupled 
time-dependent oscillators 

^oW=Eka + ^^aW0L {U^l{t)=k' + m\t)). (8) 
aa 

Then quantum states of the pion field are found by the time-dependent creation and anni- 
hilation operators [7, 8] 

aaait) = i[flait)^aa - 'fi*aait)^aa], (9) 

where itaa and (paa are Schrodinger operators. Note that these operators do not diagonalize 
the Hamiltonian (8), but satisfy the Liouville-von Neumann equations 

^Ol^lait) + [aUt), Haait)] = 0, Z-a„<,(t) + [aaait), Haa{t)] = 0, (10) 

which lead to mean-field equations for the auxiliary field variables 

^aa{t)+col{t)ipaa{t)^0. (11) 

In fact, these operators satisfy the equal time commutation relations 

[aaa{t): alp{t)] = SabSa/3: (12) 

with the aid of the Wronskian condition 

'fla^aa " ^laV^aa = i- (13) 

The Fock space for each pion field mode consists of number states defined as 

Naa{t)\naa,t)o = aL(i)«aa(^) l^^aa, ^)o = ^^aa l^aa, t)o- (14) 

We should note that these are exact quantum states of the time-dependent Schrodinger 
equation for the Hamiltonian (6) or (8). The quantum state of the pion field itself is then 



a product of each mode state for each pion field. Of a particular interest is the Gaussian 
vacuum state of the pion field 

|0,t)o = niO-'^)o- (15) 

aa 

Now the Green function for the pion field is simply given by 

Go(x,t;x',t') = nGoaa(0aa,t;0la,O, (16) 
aa 

where the (aa)-mode Green function takes the form 

Goaa{4>aa-,i'i'P'aai'^') = '^X4> aa I '^aa ) t)ooKa,t'|0lJ. (17) 

To study formation of domains during a nonequilibrium quench process, we consider the 
smooth finite quench (3). The free field theory (6) is then exactly solvable [7], and is a 
good approximation for the linear sigma model as long as |m^(t)| is larger than |0(t)| and 
(tt)^ = Y,h(} ^b/s^bp- Far before the phase transition, each mode is stable and oscillates around 
the true vacuum with the solution 

if^^iit) = ^^e-*"-*, cuai = Vk^ + mf. (18) 

V^^ai 

The two-point correlation function for each component of pion field is the Green function 
at equal times 

G'oa(x, x', t) = (0„(x, i)0„(x', t))o = G'oa(x, t; x', t) (19) 

with respect to the Gaussian vacuum or thermal equilibrium. 

On the other hand, after the phase transition (m^ = —m'j = —{rriQ — ml)), the long 
wavelength modes with k < ruf become unstable and exponentially grow, whereas short 
wavelength modes with k > ruf are still stable and oscillate around the false vacuum. Far 
after the phase transition, the unstable long wavelength modes have the asymptotic solutions 

Here /i^ and i/k depend on the quench process and satisfy the relation |//kp — |i^kp = 1- The 
correlation function is then dominated by the exponentially growing part 

Goa/(x,t;x',t) / 7^e^"-(— ')|/.kr:^^^ ^. (21) 



(27r)3 2{k^-m})' 
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Using the exact solutions [7], we obtain the two-point thermal correlation function for 
each pion field at the intermediate stage of the quench (— r < i < r) 

GoaTir,t) ~ GoaT{0,t) ^ exp(-^). (22) 

mo ^ mo 

It follows that domains obey a scaling relation for the correlation length 

Ut)-2{^y". (23) 

The power 1/4 has been found in numerical simulations [17]. At the later stage far after 
the quench (i » r) , domains still show the Cahn- Allen scaling behavior but with a different 
power 

Ut)='2{^y^\ i=t-^m)-l]{m^ + m}). (24) 

The scaling relation for the instantaneous quench is obtained by letting r = in Eq. (24). 
A kind of resonance has also been observed in the correlation function with simple poles at 

^-^2 (n = l,2,3,---). (25) 



[m 



k2)l/2' 



This structure implies certain adjusted quench rates r may lead to sufficiently large domains 
[7]. 



IV. SCATTERING, DECOHERENCE AND ORDER PARAMETER 

The tree level approximation in Sec. 3 does take into account neither any interaction 
among different modes of each pion field nor the interaction between pions. Similarly the 
Hartree approximation includes only mean-field effects among modes and pion fields and 
thus neglects any direct scatterings among modes. To go beyond the Hartree approximation 
we may use the formalism in Ref. 8. The wave functional for the Schrodinger equation can 
be expressed in terms of the Green function (kernel or propagator) as 

*(x, t)^ J ^(x, t; xo, to)*(xo, to)dxodto. (26) 

As the linear sigma model is nonlinear, we use a perturbation method. We divide the 
Hamiltonian 

H(t)^Ho(t)+XHp(t), (27) 
7 



into an exactly solvable (quadratic) Gaussian and a perturbation part 

Hp = \{r<t>af - \{4^'h)r<t>a. (28) 

Then in terms of the Green function for Hq in Sec. 3, 

- ^o(x, t))Go(x, t- x', t') = 5(x - x')5(t - t'), (29) 

we write the wave functional as 

*(x, t) = ^o(x, t) + xj Go(x, t; x', t')Hp{x', t')*(x', t')dx.'dt', (30) 

and finally obtain the wave functional of the form 

*(1) = *o(l) + A J 6-0(1, 2)^p(2)*o(2) 

J J Go(l, 2)^p(2)G'o(2, 3)^p(3)*o(3) + • • • , (31) 

where (i) denotes (xi,tj) and is a wave functional for ^o- This method goes beyond the 
Hartree approximation and gives us the wave functional in a series of A 

*(x,t)=^o(x,t) + E-^"*o"^- (32) 

n=l 

The term ^g"' comes from nth order correction of Hp. 

Wc now give a few remarks on the effects on DCC formation of the non-Gaussian (be- 
yond the Hartree) approximation. First, the nonlinear correction due to ^g"'' enhances the 
correlation length for domains by a factor [9] 

^ = (2n+ 1)1/1 (33) 

Higher order correction terms begin to grow provided that the duration of instability is long 
enough before crossing the inflection point. This condition may be provided by the quench 
time scale that is comparatively large but still smaller than relaxation time scale. Under 
this condition, significantly large DCC domains may lead to observable effects in heavy ion 
collisions or high energy cosmic rays. Second, the higher order correction terms in Eq. (32) 
include direct scattering effects. The direct scattering can be shown obviously in Eq. (31), 
where the Green function is simply given by 

oo 

Go(x, t; x', f) = ^ ^,(x, t)%{^, t') (34) 
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with being the Fock states of Hq. Thus the nonhnear perturbation Hp scatters "^aaq 
into '^b/3p, and vice versa. In particular, direct scatterings with short wavelength modes 
(environment or noise) would lead to decoherence of long wavelength modes. Therefore, 
long wavelength modes achieve not only classical correlation but also decoherence [18] and 
long wavelength modes emerge as a classical order parameter. 

V. CONCLUSION 

The QCD with two massless quarks would undergo a second order phase transition. 
The effective theory for the QCD phase transition is the hnear sigma model for the quark- 
antiquark condensate, that is, the sigma and pion fields. Under a rapid cooling process, 
the QCD phase transition proceeds far from equilibrium (nonequilibrium). In this talk we 
focused on the dynamical process of the nonequilibrium phase transition and its implications 
on DCC domain growth. 

The most prominent feature of nonequilibrium second order phase transitions is the in- 
stabihty of long wavelength modes. These modes begin to grow exponentially at the onset 
of phase transition while rolling over the barrier from the false vacuum to the true vacuum. 
Therefore, this nonequilibrium dynamical process leads to large domains of disoriented chiral 
condensate, regions of misaligned vacuum. On the contrary, domains formed from thermal 
equilibrium have small sizes determined by thermal energy and are randomly oriented in 
isospin space. 

Using the nonequilibrium linear sigma model with a smooth finite, we observed that 
DCC domains grow according to some power-law scaling relations. We found the power 

^1/4 

for the scaling behavior in the intermediate stage and the Cahn-Allen scaling power 

tV2 at 

the later stage of phase transition. However, in the Hartree approximation, exponentially 
growing long wavelength modes contribute A(7r^), which in turn attenuates the instability. 
The instabihty completely stops when X{n'^) dominates over mj. Therefore, domains grow 
for a limited time and reach a typical size of 1 ~ 2 fm in the Hartree approximation [14]. 
This size of DCC domain may not be large enough to lead to any significant observation. 

It was shown that higher order quantum corrections increase the domains size by addi- 
tional factors (2n -|- 1)^/^ if the duration of instability is long enough to make higher order 
terms grow. The longer is the quench time, the larger domains are. However, relaxation 
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time scale gives a limitation on domain growth through instability, because thermal equili- 
bration competes with instability when the quench time is comparable to relaxation time. 
This non-Gaussian effects on DCC domains may lead to observations in heavy ion coUisions 
and high energy cosmic rays. This mechanism should be distinguished from the anomaly 
enhanced domain formation [19]. 
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